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Abstract 

We consider an ideal gas of massive hadrons in thermal and chemical equi- 



("^ I librium. The gas expands longitudinally in an adiabatic way. This evolution 

for a baryonless gas reduces to a hydrodynamic expansion. Co 
is parametrized by the sound velocity. The sound velocity is 
dependent and is strongly influenced by hadron mass spectrum. 



'^ ' for a baryonless gas reduces to a hydrodynamic expansion. Cooling process 

is parametrized by the sound velocity. The sound velocity is temperature 
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1 Introduction 

In our last paper [1], there have been presented results for J/'^ suppression in an 
equilibrium hadron gas. The gas expands longitudinally according to the Bjorken 
pattern [2] and its initial conditions have been estimated from the data [3,4]. Since 
our purpose in Ref.l. was to study J/^ suppression patterns, we left aside more 
detailed description of properties of the hadron gas subjected to the longitudinal 
cooling. Now, we would like to discuss this subject more carefully. 

In the case of hadron-hadron collisions we shall assume an existence of a " central 
region" in the rapidity variable [2]. We are leaving aside here a problem of an 
appearance of a quark-gluon plasma phase. We begin our analysis with the formation 
of an equilibrium hadron gas, subjected to the longitudinal expansion. This gas 
consists of different species of hadrons and it has a non-vanishing baryon-number. 
This fact corresponds to present experimental conditions. A neutral central region 
is expected for higher then CERN collision energies of heavy ions. 

2 Equation of state 

An equation of state is given by an assumption of a thermal and chemical equilibrium 
of an ideal hadronic gas. Particle ratios are given by the temperature and chemical 
potentials related to conserved quantum numbers - strangeness and baryon-number. 
In our model we take into account all species of hadrons up to Vt^ baryon. 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of / 
species of particles, energy density e, baryon-number density n^, strangeness density 
ns and entropy density s read {% = c= 1 always) 

^^ i=i -^0 eicp^-hf^^ +gi 

^^ i=i -^0 exp^-Y^j+gi 



where Ei = [mf +p^)^/^ and mj, Bi, Si, jj,i, Si and gi are the mass, baryon-number, 
strangeness, chemical potential, spin and a statistical factor of specie i respectively 
(we treat an antiparticle as a different specie). 

And fii = BifiB + Sifis, where /i^ and ns are overall baryon-number and strange- 
ness chemical potentials respectively. 

These equations can be supplemented by a similar equation related to the electric 
charge density. A chemical potential ^q would be introduced due to the electric 
charge conservation. However, it is easy to see that for the ratio A/Z = 2 this 
chemical potential is always equal to zero. So we can omit an equation for the 
electric charge density at present experiments. 

3 Time evolution 

In general, all densities of the left sides of eqs.la-d are functions of time during the 
cooling of the hadron gas. This means that gas parameters such as temperature 
and chemical potentials should be functions of time. We would like to obtain these 
functions explicite. This can be done by solving numerically the system of equations 
for entropy density s, baryon-number density Ub and strangeness density ns, with s, 
Hb and ns given as time dependent quantities. Strangeness density is equal to zero 
during all the evolution process, but for other quantities some dynamical assumption 
is needed. In the following we shall assume a longitudinal expansion and we shall 
neglect transversal degrees of freedom. 

From the baryon-number conservation we have then 

UBit) = ^ , (2a) 

From the Bjorken model with vanishing baryon-number density we have the follow- 
ing solution for the longitudinal expansion [1,5] 

sit) = '-f (26) 

where Sq and n^ are initial densities of the entropy and the baryon-number respec- 
tively. This gives an entropy conservation during the expansion stage. We shall 
assume an adiabatic cooling process, so eq.2a will be taken for granted also with 
nonvanishing baryon-number density. This means a deviation from a simple hydro- 
dynamical model. 



To solve (lb,c,d) with s and hb given by (2) and ns = 0, we need to know 
initial values sq and n^. To estimate initial baryon- number density n% we use 
experimental results of [3]. These results are for S-S collisions, but since there are 
no data on baryon multiplicities for heavier nuclei we have to evaluate them in some 
way. 

We assume that the baryon multiplicity per unit rapidity in the CRR is propor- 
tional to the number of participating nucleons. For a sulphur-sulphur collision we 
have dNs/dy = 6 [3] and 64 participating nucleons. For an 0-U collision we can 
roughly estimate the number of participating nucleons at 16 + 58 = 74. 

The second factor of the sum has been obtained by the following assumption: 
since an oxygen nucleon is much smaller than an uranium one, we can approximate 
the part of the uranium, through which the oxygen passes, by the cylinder of the 
volume equal to ttRq ■ 2Ru. The same procedure can be apphed to the S-U case. 
Here we obtain 32 + 93 = 125 participants. 

Therefore, we have dNs/dy = 7 and dNs/dy = 11.7 for 0-U and S-U collisions 
respectively. Having taken the initial volume in the CRR equal to 7tR\ • 1 fm, we 
arrive at n^ = 0.25 fm~^ for both cases. 

To find Sq) first we have to solve (la,b,c) with respect to T, /i^ and /i^, where we 
put e = Co and so on. For cq we have taken estimates given in [4]. As a result, we have 
obtained Tq = 212 MeV, ^i% = 34.6 MeV and 12% = 133 MeV for the 0-U colhsion 
(eo = 2.5 GeV/fm^) and Tq ^ 209.1 MeV, /x^ = 36.7 MeV and fi% ^ 143.3 MeV for 
the S-U collision (eg = 2.3 GeV/fm^). Then, from (Id) we have Sq = 13.68 fm^^ for 
0-U and sq = 12.74 fm'^ for S-U. 

Now, having put (2) and ns = into (lb,c,d), we can solve them numerically to 
obtain T, fis and /x^ as functions of time. 

4 Results 

Our results are presented in fig.l, 2, where solid, long-dashed and dashed lines 
mean the temperature, the strangeness chemical potential and the baryon chemical 
potential respectively. Fig.l shows results for 0-U collision initial conditions and 
fig. 2 for S-U ones. The time scale is chosen in a way which enables the temperature to 
reach the freez-out at 140 MeV (figs. la and 2) or 100 MeV (fig.lb). This corresponds 
to the freez-out time equal to tj.o. = 10.4 fm or tfo. — 58.6 fm respectively. The most 
interesting feature of our results is the behaviour of the temperature. For figs, la. 



lb, and 2 the following temperature approximations hold respectively: 

/1\ 0-178 

T(t) ^212.4- f-j (3a) 

/I ^ 0.189 

T(t) ^217.2- (-) (35) 

/ 1 ^ 0.179 

T(t) ^209.7- (-] . (3c) 

We can see that all above expressions have the form known from the solution for the 
longitudinal expansion of a baryonless gas with the sound velocity constant, namely 

[5] 

nt) = To • (-) \ (4) 

where Cg is the sound velocity and we put the initial time to equal to 1 fm. 

We have checked that for hb = results for the temperature function are very 
similar to those in eq.3. This is shown in figs. 3 and 4. Fig. 3 has the same initial 
conditions as fig. lb and fig. 4 the same as fig.2. The following approximations of the 
temperature function hold 

/I ^ 0.180 

T(t) =215.2- (-j (5a) 

T(t) =208.9- f-j . (56) 

We can see that the power in eqs.3 and 4 is around 0.18. Therefore a question 
arises: is this value connected in any way with the sound velocity of the hadron gas? 
For the baryonless case the answer is positive. We have checked this computing 
straightforward 

which is the value of the sound velocity squared for a baryonless gas [5]. The 
results are presented in fig.5. For comparison, the square of the sound velocity of 
a pion gas is also depicted (dashed lines). We can see that in the range of the 
temperature 200-140 MeV the square of the speed of sound equals 0.17-0.18 indeed. 
We have also found a region of the temperature where the sound velocity decreases 
as the temperature increases. 

Nevertheless, the condition for the stability of the expansion [6], d/dT{sCs/T) > 
is still valid because sCs/T is an increasing function of temperature in the whole 
temperature region. This function is presented in fig.6. 
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Figure Captions 

Fig. la. Dependence of the temperature (solid line), the strangeness chemical po- 
tential (long-dashed line) and the baryon chemical potential (dashed line) on 
time. Initial conditions are chosen for 0-U collision with freez-out at 140MeV. 

Fig. lb. Same as fig. la but for freez-out at 100 MeV. 

Fig. 2. Same as fig. la but for S-U collision. 

Fig. 3. Dependence of the temperature on time for hadron gas with non-vanishing 
(solid line) and vanishing (dashed line) baryon number. Initial conditions are 
chosen for 0-U collision with freez-out at lOOMeV. 

Fig. 4. Dependence of the temperature on time for hadron gas with non-vanishing 
(solid line) and vanishing (dashed line) baryon number. Initial conditions are 
chosen for S-U collision with freez-out at 140MeV. 

Fig. 5. Baryonless gas: dependence of the sound velocity squared on time. Solid 
line: all hadrons up to Q~ resonance. Dashed line: pion gas. 

Fig. 6. Stability factor sCg/T for a baryonless hadronic gas. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411250v2 
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